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The Dynamical State of the Starless Dense Core FeSt 1-457: A 

Pulsating Globule? 

E. D. Aguti 1 ' 2 , C. J. Lada 2 , E. A. Bergin 3 , J. F Alves 4 , and M. Birkinshaw 1 

ABSTRACT 

High resolution molecular line observations of CS (J = 2 — > 1), 
HCO+ (J = 1 — ► 0), C 18 (J = 1 — ► 0), C 18 (J = 2 -f 1) and 

N2H + ( J = 1 — > 0) were obtained toward the starless globule FeSt 1-457 in 
order to investigate its kinematics and chemistry. The HCO + and CS spectra 
show clear self-reversed and asymmetric profiles across the face of the globule. 
The sense of the observed asymmetry is indicative of the global presence of 
expansion motions in the outer layers of the globule. These motions appear to be 
subsonic and significantly below the escape velocity of the globule. Comparison 
of our observations with near-infrared extinction data indicate that the globule 
is gravitationally bound. Taken together these considerations lead us to suggest 
that the observed expansion has its origin in an oscillatory motion of the outer 
layers of the globule which itself is likely in a quasi-stable state near hydrostatic 
equilibrium. Analysis of the observed linewidths of C 18 and N 2 H + ( J = 1 — > 0) 
confirm that thermal pressure is the dominant component of the cloud's internal 
support. A simple calculation suggests that the dominant mode of pulsation 
would be an 1 = 2 mode with a period of ~ 3 x 10 5 yrs. Deformation of the 
globule due to the large amplitude 1 = 2 oscillation may be responsible for the 
double-peaked structure of the core detected in high resolution extinction maps. 
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Detailed comparison of the molecular-line observations and extinction data 
provides evidence for significant depletion of C 18 and perhaps HCO + while 
N 2 H + ( J = 1 — > 0) may be undepleted to a cloud depth of ~ 40 magnitudes of 
visual extinction. 

Subject headings: ISMxlouds, ISM: globules, ISM: molecules, ISM: individ- 
ualist 1-457) 



1. Introduction 



A key component to understanding star formation lies in determining the dynamical 
nature of globules and cores of molecular clouds from which stars form. This is typically 
done using molecular lines, and observations of asymmetry or self-reversal in the line profiles 
of optically thick species portray the motion of the outer layers of material within a cloud. 
Many of the spectra measured from starless cores and globules have blue asymmetirc profiles, 
such that the blue side of the profile is stronger than the red side, suggestive of r ed-shifted 



self-a bsorption and infall motion. For example, in their survey of starless cores, iLee et al. 



( 120011 ) found that approximately one third of the cores they observed displayed bona fide 
infall motions in the outer layers and therefore were consistent with expectations of collapse. 
In this survey and in a number of others, however, there were a number of cores that displayed 
evidence for blue-shifted self-absorption in their spectra, e.g., CB246, L1521F and L429-1. 
This is an indication of expansion but with no protostar in the core, what mechanism drives 
the expanding motions? There are also cores and globules which show both redshifted and 
blueshifted self-absorption simult aneously in a s ingle map of one molecular transition, for 
example, L1512, L1689B, L1544 Lee et aliboOlL and B68 Lada et al.l l2003h . Other cores 
and globules have spectra that shift from redshifted to blueshifted self-absorption depending 
on which molecular transition is observed. For example, self -absorption in the source LI 157 
is redshifted in the HCO+fJ = 3 - > 2) line profile but blueshifted in the HCO+(J = 4^3) 
line profile (IGregersen et al.lll997l ). Clearly complex motions are occurring in the various 
layers of the cores but it is not clear what causes them. 



Redman et al.l (120041 ) concluded that the pattern of redshifted and blueshifted self- 



absorption observed across the source in L1689 is caused by rotational motion in the globule. 
In B68, where the spatial pattern switches from blueshifted to redshifted profiles both east 
and west of the cent e r of t he globule, the observed line profiles cannot be explained by 
rotation. lLada et al.l (120031 ) proposed that B68 could be in a quasi-stable state near hy- 



drostatic equilibrium, and the systematic velocity field observed is likely to be caused by 
oscillatory motions within the outer layers of the globule. For stable and quasi-stable clouds, 
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such motions are predicted by theory, e.g.. iMatsumoto h Hanawal (120031 ) showed that small 
amounts of rotation can cause a n initially gravitationally unstable cloud to stabilize and 
oscillate; Oarke &: Pringld (119971 ) in considering the thermal and dynamical balance in low 
mass dense cores find that a quasi-equilibrium state can be established, whic h is not static 
but is pulsationally stable, and in their hydrodynamic study of starless cores, iKeto fc Field 
( 120051 ) found that cores could oscillate with periods of about one million years (or the sound 
crossing time) if perturbed by a modest external force. Such oscillatory behavior is a signa- 
ture of dynamical equilibrium, and therefore, the presence of such complex patterns in their 
spectral line maps could indicate pulsationally stable initial states for the star-forming cores 
and globules. 

In this stud y, we have obtained high angular resolution spectra of the source listed in 
the catalogue of iFeitzinger fc Stuwd (119841 ) as FeSt 1-457. It is a small, round, dark globule 
situated in the Pipe Nebula with no associated IRAS, MSX or Spitzer Space Telescope 
point sources, and is therefore almost certainly starless. Its distance has not been directly 
determined but it is associated with the Pipe Nebula, for which some distance estimates 
have been made. The most reliable estimate t o date is 130 pc deter mined by using infrared 
extinction measurements for Hipp arcos stars, (ILombardi et al.ll2006T) . In their Bonnor-Ebert 
sphere fitting analysis, however, iKandori et al. (boosh report a distance of 70 pc, which 
compares to the value of 90 pc derived in a similar fashi on for B68, also, incide nt ly, part of 
the Pipe Nebula complex. In this paper, we adopt the ILombardi et al.l (120061 ) distance of 
130 pc. The angular diameter of FeSt 1-457, as apparent in optical images, is ~ 4.5', which 
is ~ 35,000 AU (0.17 pc) at 130 pc. Section [2] describes our observations of FeSt 1-457, and 
^31 is a description of the results. We detected expanding motions in HCO + and CS line 
profiles across the surface of FeSt 1-457 but since no protostar has been detected, the source 
of this expansion is a mystery. In §H we show that the core is bound and self-gravitating 
and we propose that pulsation or oscillatory motion within the cloud layers may account for 
the observed expanding motions. Finally, we present evidence for molecular depletion in the 
central regions of the core. 



2. Observations 

The observations reported here were made in July 2003 using the 30-m IRAM millimeter- 
wave telescope at Pico Veleta in Spain. The dual mixer, dual-channel receiver was tuned to 
observe the molecules listed with their respective frequencies, beamwidths and velocity reso- 
lutions in Table [TJ The frequency- switching mode was used to obtain the observations, and 
system temperatures were 150 - 200 K. The chopper wheel method was used for calibration. 
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The spectrometer is an autocorrelator configured to give the velocity resolutions listed in 
Table[TJ Beam efficiencies were ~ 0.75 for the lower frequency 80 - 115 GHz range, and ~ 
0.6 for the higher frequency 197 - 266 GHz range. Observations were made in a grid with 
a reference position at a 200 o = 17 /l 35 m 47!5, 5 200 o = -25°33'2?0. HCO+, C 18 ( J = 1 -> 0) 
and C 18 (J = 2 — > 1) were uniformly sampled such that a region centered on the reference 
position was observed every 24" out to a maximum distance of 96 " in Right Ascension and 
120" in Declination, a total of 9 x 11 or 99 spectra for each species. This grid spacing is ap- 
proximately one beamwidth, slightly over for HCO + and under for C 18 0. The N 2 H + and CS 
emission was much weaker and less extended so the grid size was reduced to 5 x 5 around 
the reference point with an additional two measurements directly north, east, south and 
west of the center, resulting in a total of 33 spectra for these species. The spacing between 
measurements was also 24", approximately one beamwidth for each of the molecules. Data 
reduction was done using the standard IRAM package CLASS and IDL (Interactive Data 
Language by Research Systems, Inc.). 

3. Results 

3.1. Self-Reversed HCO+ and CS Lines 

Figured] shows the CS (J = 2 1), HCO+ (J = 1 -> 0) and C 18 ( J = 1 -> 0) spectra 
at the reference position. The CS and HCO + profiles display a split asymmetric, double- 
peaked shape while the C 18 line is single-peaked, though slightly asymmetric. The C 18 
line appears to bisect the CS and HCO + profiles in velocity indicating that the two latter lines 
are likely optically thick and self-reversed. The blue-shifted sense of the self-absorption in the 
self-reversals is indicative of expansion motions in the outer cloud layers. Individual spectra 
of the HCO + and C 18 ( J — 1 — > 0) emission from across FeSt 1-457 are simultaneously 
displayed in Figure [2] for comparison. Examination of the figure shows similar double-peaked, 
self-reversed HCO + emission profiles at several positions around the center of the globule. 
In all cases the sense of the asymmetry is the same suggesting expansion of the outer cloud 
layers. In most other positions the HCO + lines display blue-shifted self-absorption relative 
to C 18 0, also consistent with expansion motions. 

As can be seen in Figure |3]CS spectra are also clearly self- reversed with blue-shifted self- 
absorption similar to the HCO + lines in positions where the HCO + was also self-reversed, 
but the signals have lower intensity. Similar to HCO + , the other CS spectra appear to 
exhibit asymmetric profiles relative to C 18 with a sense indicative of expansion motions. 
The critical density (3 x 10 5 cm -3 ) at which CS emission occurs is the same as that for 
HCO + emission (Ungerechts et al. 1997) so it is reasonable to assume that the two emission 
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lines probe the same layer of material in the cloud. 

A rough estimate of the expansion speed of material was determined by comparing 
the velocity of the stronger HCO + peak to the peak velocity of a C 18 spectrum at the 
same position. The peak velocities were obtained by fitting Gaussian line profiles to the 
spectral lines. This resulted in velocity differences SV = v p k{C 18 0) - v p k(RCO + ) which are 
all negative confirming that the blueshifted profiles are characteristic of expansion, and 
indicating that material is expanding in the cloud over the whole layer under observation 
with a mean expansion velocity of -0.09 ± 0.04 km s" 1 . The same process was applied to 
the CS spectra and the velocity differences (SV = v p k(C ls O) - v p k(CS)) were also found to 
be negative everywhere with a mean difference (expansion velocity) of -0.12 ± 0.02 km s^ 1 . 
This is similar to the range of SV for HCO + . This expanding motion is therefore also evident 
in the CS line emission profiles. 

Another estimate of the exp ansion speed of material within the cloud was obtained using 
the model of iMyers et al.l (119961 ). This model can only be applied to positions (10 positions) 
where there is a clear double peak. In this model, the expansion velocity, v exp is given by, 

a 2 1 + e(T BD /T D ) 

V exp ln frr Irri \1 V / 

vr-vb 1 + e(T RD /T D ) 

where Tp is the brightness temperature of the dip, Tbd is the height of the blue peak above 
the dip with its corresponding velocity, vb, T rd is the height of the red peak above the dip 
with its corresponding velocity, Vr, a is the velocity dispersion of an optically thin line (C 18 
here). For the HCO + lines, the mean expansion speed was calculated to be -0.07 ± 0.02 and 
and for the CS lines, -0.07 ± 0.02 km s~ l ; both these estimates are somewhat lower than those 
derived from the peak velocity difference method in the previous paragraph. Nonetheless, 
though somewhat uncertain, the overall expansion speeds we estimate are formally less than 
the one dimensional sound speed of a ~ 0.19 km s _1 in a 10 K gas. 



3.2. Velocity Field 

Figure H] shows maps of the variation of the velocity of the peak in the C 18 spectra. 
The peak velocities were determined from Gaussian fits to the line profiles. We note that 
a few of the C 18 ( J — 1 — > 0) and C 18 (J = 2 —>■ 1) lines are flat-topped, broadened or 
slightly asymmetric, indicating that the lines are slightly optically thick. The two maps, 
however, do reveal a systematic velocity gradient increasing from upper left to lower right of 
the map. 



To estimate the magnitude and direction of this gradient, the method of 



Goodman et al 
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( 119931 ) was used as adapted by lLada et al.l (120031 ). The velocity gradient is assumed to be 



linear when projected against the plane of the sky so the observed velocity vi sr can be related 
to the velocity gradient dv/ds using 

dv dv 
visr = vo + —Aa cos9 + —A5 sm9 (2) 
ds ds 

Aa and AS are Right Ascension and Declination offsets in arcseconds. vq is the systemic 
velocity of the cloud and 9 is the angle between north and the direction of the velocity 
gradient of magnitude dv/ds. A least-squares fit of a two-dimensional plane to the ob- 
served Vi sr (line center velocity) measurements of C 18 ( J = 1 — > 0), C 18 ( J — 2 — > 1) and 
N2H + ( J — 1 — > 0) (the N2H + ( J — 1 — > 0) spectra contained hyperfine components, and so 
the isolated component at the same vi sr as the C 18 ( J = 1 — > 0) line was used in this anal- 
ysis) provided estimates given in Table [2j The errors for the N 2 H + are larger because there 
were only nine spectra measured that were useful for this purpose. The N 2 H + emission 
was less extended than the C 18 emission, therefore, more densely sampled observations are 
needed to confirm the estimates using the N2H + data. 

If the velocity gradient derived previously is removed, the average velocity difference 
between the neighboring spectra is then essentially zero (~ 0.0025 km s _1 ). 



3.3. Line Width Distribution 

Line widths (FWHP) were determined from Gaussian fits to the observed C 18 (J = 1 — > 0) 
and C 18 (J = 2 — > 1) line profiles. The N2H + ( J = 1 — > 0) spectra have several hyperfine 
components so the line widths were determined from a Gaussian fit to the isolated (Fi = 0—1) 
component. The resulting C 18 ( J = 1 — ► 0) and C 18 ( J = 2 — > 1) line widths were found 
to be in a range from ~ 0.19 to ~ 0.35 km s _1 . The N2H + line widths were narrower, and 
were in the range ~ 0.15 to ~ 0.25 km s" 1 . Figure [5] shows the variation in line width of the 
C 18 ( J = 1 — > 0) line profiles. Because some of these lines are optically thick (see comment 
in § 13.21) . the line widths determined from a single Gaussian fit are slightly larger than they 
would be if the lines were optically thin. Nevertheless, the line widths seem to increase 
slightly at positions away from the highly extincted region around offset (24", 0"), marked 
with a white cross in Figure [5j This is similar to B68 where the linewidths also appeared to 
increase away from its center. The reason for this is not clear. 

No independent measure of the kinetic temperature of FeSt 1-457 has as yet been made; 
a value of 10 K has therefore been adopted in th is paper because data has shown that 



this applies to most globules and dense cores (e.g., iBenson fc Myersl (119891 ) ) . The thermal 



line width, v /[(81n2)/cT K -/(/xmj/)] for C 18 ( J — 1 — > 0) lines was therefore calculated to be 
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0.123 km s 1 , and for N 2 H + lines, 0.125 km s 1 . The nonthermal contribution to the line 
width was then estimated using the relation, 

(Av obs ) 2 = (Av th ) 2 + (Av nth ) 2 , (3) 

where Av b s is the observed line width, Av th is the thermal line width and Av nt h is the non- 
thermal line width. The resulting average nonthermal line width for the C 18 ( J = 1 — > 0) 
lines was ~ 0.25 km s" 1 , and for the N 2 H + lines, ~ 0.15 km s _1 . To the extent that these 
lines are optically thin, this suggests that in producing the observed line profiles, nonthermal 
broadening mechanisms, e.g., turbulence, may play a role. There may be more turbulent 
motion in the outer layers probed by the C 18 molecular lines than in the inner more dense 
layers probed by the N 2 H + line. However, the corresponding one dimensional non-thermal 
velocity dispersions (<T nt h) are 0.11 and 0.06 km s _1 for the C 18 and N 2 H + emitting gas, 
respectively. These values are both subsonic, that is, less than the one dimensional sound 
speed (0.19 km s -1 ) in a 10 K gas. 

3.4. Spatial Distributions of Molecular Emission and Dust Extinction 

In Figure 6 we show contour maps of C 18 (1-0), N 2 H+ (1-0), and HCO+ (1-0) inte- 
grated emission overlaid on a grey-scale map of the distribution of visual extinction. The 
extinction map was constructed from the data of Alves et al. (2002) and convolved with 
a 24 arc sec Gaussian smoothing kernel to match the resolution of the molecular-line data. 
The patterns seen in this figure are comparable to those seen in numerous other low-mass 
pre-stellar cores such as B68 or L1544 (Bergin et al. 2002; Tafalla et al. 2002). In particular, 
both C 18 and HCO + show emission distributions that are broader than the distribution in 
extinction with peaks well separated from the extinction maximum. In contrast N 2 H + shows 
the highest degree of agreement with the dust extinction. This pattern is attributed to the 
freeze-out of molecules on the surfaces of cold dust grains in gas where the density exceeds 
10 5 cm" 3 (e.g. Bergin & Langer 1997; Aikawa et al. 2005). In this scenario the freeze-out 
of CO and its isotopologues leads to the formation of N 2 H + in the dense gas dominated by 
freeze-out and traced by the dust. HCO + has structure similar to CO which is not surpris- 
ing since HCO + likely forms from gas phase reactions involving CO. For a summary of this 
process in starless cores see Bergin & Tafalla (2007). 
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4. Discussion 

4.1. Is FeSt 1-457 Gravitationally Bound? 

One explanation for the expansion of material from FeSt 1-457 could be that the globule 
is simply unbo und. This might be expected if the core is a transitory feature in a global 



turbulent sea (IBallesteros-Paredes et al.ll2006l ). In this section, we assess this possibility by 



determining whether or not the core is bound. A mass for FeSt 1-457 of ~ 3.1 [^] 2 M was 
derived by spatially integrating the extinction data over the area of the globule, after cor- 
recting for a significant foreground/background extinction component (Ay ~ 6 magnitudes) 
and an assumed distance of d pc. The magnitude of the background extinction was derived 
from examination of the azmuthially averaged extinction profile of the globule constructed 
from the Alves et al. (2002) deep extinction data and is the same as that derived by Alves 
et al. (2007) for this core from a wavelet decomposition of their 2MASS extinction map of 
the Pipe cloud. The escape velocity (y/[2GM/R\) from FeSt 1-457 is estimated to be ~ 0.6 
[^io ]°' 5 k m s_1 - The average three-dimensional velocity dispersion (a/3<2 2 + 3cr^ th ) of the 
bulk gas in the globule is ~ 0.3-0.4 km s _1 , and significantly less than the escape velocity. 
Thus the globule is likely to be gravitationally bound for the assumed distance of 130 pc 
or for any other plausible distance to the cloud. Moreover, the expansion speeds detected 
in the self-absorbed molecular line data (~ 0.1 km s -1 , see § 13. ip are also significantly less 
than the calculated escape speed. The expansion of the outer layers of the globlule appears 
not to be caused simply by evaporation of its gas as in an unbound, transient cloud. 

A Jeans mass of 3.8 M was derived using Mj = 18 M Tr 1 ' 5 n -0,5 where T is 
the globule's kinetic temperature assumed to be ~ 10 K. However, if we assume the lines 
are only thermally broadened, with no turbulence, then the kinetic temperature is 17 K 
and this doubles the Jeans mass estimate (§ 13.31) . The average density n=2.2 x 10 4 cm -3 
was determined using the extinction data and a radius, R = 2.6 x 10 17 cm. The radius 
was determined to be the point at which the radially averaged extinction profile reaches 
the background and marks the outer edge of the globule. Since the mass of FeSt 1-457 
is comparable to its Jeans mass, we conclude that it is likely marginally stable against 
gravitational collapse. If thermal pressure is not the only source of internal support then the 
globule could be even more stable against collapse. More detailed analysis of the globule's 
structure would be required to better evaluate its overall stability. 
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4.2. Possible Rotation? 



There may be some contribution to stability of the globule from rotation. Assuming 
soli d body rotation, 0, th e ratio of rotational kinetic energy to gravitational energy, is given 
by (IGoodman et al.lll993l ). 



P 



{l/2)Iu 2 _ lpu 2 R i 
qGM 2 /R ~ 2q GM 



(4) 



R is the globule's radius, M its mass and p is defined such that the moment of inertia, 
I = pMR 2 , and q is defined such that the gravitational potential energy is qGM 2 /R. 
was estimated to be 0.01 for FeSt 1-457 using M = 3.1 M , R = 2.6 x 10 17 cm, uo = 
(dv/ds)/ smi = [1.5 km s _1 pc _1 ]/sini (for C 18 ( J = 1 — > 0), see Table [2]), where i is the 
inclination of the globule to the line of sight, and assuming p/q = 0.22 as for a sphere with 
a density profile r~ 2 and sinz = 1. The contribution of rotation to the overall stability of 
Fe St 1-457 is there f ore no t significant. Note that of 0.01 is also consistent with the results 
of IGoodman et al.l (119931 ) that most clouds have < 0.05. 



4.3. Thermal vs. Turbulent Support 

Thermal motions may be a significant source of pressure support for FeSt 1-457. This 
can be quantified by evaluating the ratio of thermal to nonthermal (turbulent) pressure given 
by, 

a 2 

R p = — , (5) 

a nth 

where a is the one dimensional isothermal sound speed and a nt h is the one dimensional 
nonthermal or turbulent velocity dispersion and a nt h = [At>„ i ? t / v /[81n2]]. Assuming a gas 
temperature of 10 K, the average value of R p for all the C 18 observations is 3.75 ± 1.95, 
which may be an underestimate for those lines that may be somewhat thick. The average 
value of R p for the N 2 H + data is 6.09 ± 2.07. These values indicate that the thermal pressure 
significantly exceeds the turbulent pressure in the globule, and absent strong static magnetic 
fields, is the dominant source of internal support against gravity. 

In comparing the turbulent velocity dispersion to the isothermal sound speed of 0.19 km s" 
in a 10 K molecular gas, the nonthermal motions appear to be subsonic over the surface of 
FeSt 1-457. If the N2H + observations probe a deeper layer of the globule than the C 18 lines, 
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these motions are even more subsonic in the inner layers of the globule. These considerations 
further indicate that thermal motions provide a significant source of internal pressure. 

The apparent velocity gradient in the C 18 data was calculated and removed in § 13.21 
The resulting average velocity difference between neighboring spectra was essentially zero 
(~ 0.0025 km s _1 ) so there appears to be no residual scatter in peak velocities of the spectra 
on scale sizes comparable to the beam width that could be caused by turbulent motions. 
This also is consistent with significant thermal support for the globule. 



4.4. Pulsating Globule? 

In the absence of an embedded source in FeSt 1-457 that could cause outflow of material, 
it is difficult to explain the observations reported here. In the previous section we showed 
that the globule is gravitationally bound with thermal pressure as the dominant source of 
internal support. But what causes the observed expansion motions? The facts that 1- the 
globule is bound, 2- thermally supported and 3- does not exceed the Jean's mass, suggest 
that this core is stable and not far from a state of hydrostatic equilibrium. Moreover, the 
expected (one-dimensional) virial velocity dispersion, a V i r i a i = y/[^GM/R], is 0.18 km s _1 
and is comparable to the sound speed in a 10 K gas as would be expected for a thermally 
dominated, stable core. Indeed, we derive the velocity dispersion for the H 2 gas in the core 
to be a = \J a? + a\ th ~ 0.21 km s -1 close to the predicted value. However, given that its 
outer layers are globally expanding, the globule cannot be presently in a strict equilibrium 
state. One viable explanation that can simultaneously account for all these facts, is that 
FeSt 1-457 is in dynamical oscillation around an equilibrium state. The outflowing motion 
is then part of a mode of oscillation such that its layers pulse in and out with a period of 
possibly a few hundred thousand years. We further explore this interesting possibility below. 



In their study of molecular lines observed in the globule B68, lLada et al.l (120031 ) con- 
cluded that that globule could be experiencing low order mode, nonradial oscillations of its 
outer layers. They suggest that the globule may have interacted with the shock of a super- 
nova remnant, which instigated a perturbation that set at least its outer layers in oscillatory 
motion. Figure 7 of their paper displays the real parts of the spherical harmonic functions 
for a range of low order modes for an oscillating sp here. The / = 2, m = 2 mode corresponds 



closely to the pattern observed in B68. A study by lKeto et al.l (120061 ) demonstrated that lin- 
ear perturbations on a pressure-bounded thermally supported sphere could indeed produce 
the spatial variation in profile shape observed in B68. FeSt 1-457 could also be oscillating 
in this mode but the geometry is such that we are looking 'edge-on' to the pulsating cloud. 
This means the mode of pulsation appears similar to the 'breathing' or / = mode, i.e., 



outflow motion over the whole surface of the cloud (IKeto et al 



2006|). 



Calculations were carried out to find the modes of pulsation for an isothermal, spherical 
globule of gas with similar parameters (e.g., radius = 2.6 x 10 17 cm, density = 2.2 x 10 4 
cm -3 , internal sound speed = 0.18 km s -1 ) as for FeSt 1-457. Small wave-like perturbations 
were introduced, and a differential wave equation was determined using the perturbed and 
unperturbed forms of its equations of state. The modes of oscillation postulated are likely to 
be acoustic, with relatively large amplitudes and energies, so that the linear approximation 
that we use for the mode frequencies is a simplification of the full dynamical problem. In this 
linear approximation, we ignore the gravitational modes, and find that the acoustic modes 
follow the dispersion relation 



h[ ^ (x w ) C ou t Po(out) ji {%uj [Cout/ c in]) 
Cin Po(in) jl \%ui {pout/ Qn]) 



DM 



h\ feu 



(6) 



= u)R§lc ou t where u is the frequency of the oscillations; R is the unperturbed radius; 

c ou t are isothermal sound speeds inside and outside the globule respectively; po(m) , Po(out) 
are the unperturbed densities inside and outside the globule respectively; and, and ji are 
spherical Hankel and Bessel functions of the first kind or order I, with the prime denoting 
differentiation with respect to the argument. 

The frequency uj is complex and roots of the dispersion relation where uj has a negative 
imaginary part, correspond to decaying modes of oscillation. The required mode for pulsation 
is the least-damped mode or the mode which decays at the slowest rate; this is the mode 
with the smallest negative imaginary part of uj. This mode was calculated to be I = 2 with 
a frequency of ~ 9 x 10~ 14 Hz corresponding to an oscillation period of ~ 3 x 10 5 years, 
comparable to the sound crossing time. It is plausible that this oscillation was excited 
by some transient dynamical event, such as a shock, that deposits an energy greater than 
about 3 x 10 43 ergs, (our best estimate of the kinetic energy in coherent motions) into the 
outer regions of FeSt 1-457 , where HCO + m olecules exhibit a coherent pattern of expansion 
motions. Calculations by IKeto et al.l (120061 ) show that such large amplitude motions (i.e., 
motions comparable to the sound speed) can cause significant geometrical distortions of a 
core. An im portant point to n ote is that in the high resolution extinction map of FeSt 1-457 



obtained by lAlves et al.l (120021 ). a slight splitting of the core's central regions is visible. This 
splitting could be another manifestation of an I = 2 mode of oscillation. 

FeSt 1-457 is situated in the dire ction of the Pipe N ebula, which lies at the edge of the 
Scorpio Centaurus OB Association. lOnishi et al.l (119991 ) suggested that stellar wind from 
the B0 type star r Sco caused compression of the molecular gas, triggering star formation 
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in the B59 cloud, s i tuated 
paper 



5° west and ~ 3° north of FeSt 1-457. In Figure 5 of their 
Onishi et al.l (119991 ) show an observation of FeSt 1-457, designating it 'Core 12'. In 



discussing the effect of the OB stars on the Pipe Nebula, lOnishi et al.l (119991 ) estimated that 
~ 1 x 10 46 ergs has been deposited in the cloud over a period of 1 x 10 7 years. If FeSt 1-457 
is indeed near this OB association and was also able to intercept only 0.3% of the estimated 
energy deposition, then the effects of the postulated shock waves may be more widespread 
than previously speculated. Not only did they trigger star formation in B59, but they may 
have also set gaseous globules such as FeSt 1-457 and B68 (which is also part of the Pipe 
complex) into oscillation. More observations and comparison with theoretical models of 
cloud pulsation are needed to confirm this hypothesis. 



4.5. Molecular Depletion 

In Figure 7 we provide a direct comparison of the dependence of C 18 (1-0) and N 2 H + 
(1-0) emission with visual extinction at the same angular resolution. For both species the 
figures display trends that are similar to those seen in B68 by Bergin et al. (2002). In these 
plots, a linear correlation between integrated intensity and Ay corresponds to a constant 
abundance throughout the cloud. 

The C 18 (1-0) emission shows a linear dependence with Ay until ~ 10 — 12 mag 
whereupon the trend flattens. Given the drop in the the C 18 integrated emission/ Ay 
relation near Ay ~ 12 mag we have fit the following function: J T m bdv(C l8 0) = a+b{Ay — 12) 
to the data. We find an intercept of a = 1.09 ± 0.01 K km s -1 and a slope of 

_ fo. 117 ±0.002 ifAy<12 m , 

\o.002± 0.001 ifAy>12 m . ' 

Thus for Ay < 12 m , where we see a significant linear correlation between gas emission and 
dust extinction, we estimate a C 18 abundance of ~ 10 -7 (relative to H2). This is comparable 
to that measured previously by Frerking, Langer, & Wilson (1982) of 1.7 X 10 -7 in Taurus and 
to that obtained from a similar comparison of extinction and CO measurements toward L 977 
(Alves et al. 1999). For Ay > 12 m the correlation is essentially flat, indicative of saturation 
or, as in B68, a drop in abundance towards denser regions. We can use the intercept to 
estimate average abundance in this gas and find that the abundance is ~ (2 — 4) x 10 -8 
between 20 and 40 magnitudes. Thus, we estimate that CO is depleted by a factor of ~5. In 
§3.2 we noted that the C 18 emission lines show indications of being slightly optically thick. 
Thus these abundance estimates are upper limits to the average depletion in the cloud. 
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The situation for N2H + (1-0) in Figure 7 is more complicated. There is a rough cor- 
respondence between integrated emission and Ay and it is well fit by a line with a slope of 
0.106±0.001 K km s _1 mag^ 1 and an intercept of -0.711±0.031 K km s -1 . This is consistent 
with an abundance of 10~ 10 relative to hydrogen. However, we cannot exclude the possibility 
that the trend becomes flat at Ay > 20 mag. Fits to the intensities of the hyperfine lines of 
N 2 H + (J = 1 — > 0) suggest that the total opacities are of order 12 at the emission peaks, so 
that all the hyperfine components are saturated. This is not the case at low Ay, where the 
lines appear to be optically thin. Thus we believe that the line integrated intensity-Ay cor- 
relation may be turning over at high Ay as a result of saturation and rather than depletion. 
However, we also note that the spatial map of N 2 H + in Figure 6 displays peaks northeast 
and southwest of the dust extinction peak and this could also be an indication of depletion 
of this species similar to what has been observed in B68 (Bergin et al. 2002). However, due 
to the high opacity of the line, it is not possible to separate the effects of spatial variations 
in excitation from those of column density without detailed chemical modeling. 



5. Summary 

High resolution radio observations were made of the globule FeSt 1-457. Molecular lines 
observed were CS ( J = 2 -> 1), HCO+ (J = 1 -> 0), C 18 ( J = 1 -> 0), C 18 (J = 2 -> 1) 
and N 2 H + (J = 1 — > 0). The HCO + and CS spectra showed clear self-reversed, asymmetric 
profiles across the face of the globule. The sense of the asymmetry is indicative of the presence 
of global expansion motions across the globule. These motions appear to be subsonic and 
significantly below the escape speed of the cloud. A search for IRAS, MSX sources and 
Spitzer sources showed the globule to be starless, and therefore the origins of these expansion 
motions is unclear. 

In this paper, we propose the explanation that the globule is in a quasi-stable state 
near hydrodynamic equilibrium such that its outer layers are pulsating. It is possible that 
a passing shock wave from a nearby OB star association has excited the globule setting its 
outer layers into motion. Such oscillatory motion is visible in the layers on the Sun (where 
the fundamental radial period is ~ 1 hour and other oscillation periods are ~ 5 minutes) but 
in FeSt 1-457, the perio ds are thought to be of the order of hundreds of thousands of years. 



This is similar to what lLada et al.l (j2003l ) observed in B68. Calculations to determine the 
dominant mode of pulsation resulted in an I = 2 mode with a frequency of ~ 9 x 10~ 14 Hz, 
and an oscillation period of ~ 3 x 10 5 years. A high resolution extinction map of FeSt 1-457 
exhibits a double peak structure at the h ighest levels of opacity, corresponding to visual 



extinctions of ~50 mag (lAlves et al.ll2002l ). It is possible that the proposed I = 2 mode of 
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oscillation could explain the observed splitting of the core in this column density map. 

Finally, we find evidence for significant depletion of CO and perhaps HCO + in this 
globule. However, N 2 H + ( J — 1 — > 0) may be undepleted to depths of about 40 magnitudes 
of visual extinction in the core of the cloud. 

We are grateful to Dr. Tracy Huard and the staff at the IRAM 30 m telescope for 
their assistance with the observations. We thank Dr. Carlos Roman-Zuniga for assistance 
in constructing the appropriate extinction map and with calculating the cloud mass. We 
thank the referee for insightful suggestions that strengthened the presentation. EDA is 
particularly indebted to Dr. Mike Masheder for his able guidance, useful suggestions and 
continued support throughout all aspects of this project. EDA was supported by a PPARC 
Postgraduate Studentship. CJL acknowledges support from NASA Origins grant NAG- 
13041. 
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Fig. 1.- Emission-line spectra of CS ( J = 2 -> 1), HCO+ (J = 1 -> 0) and 
C 18 ( J — 1 — > 0) near the center of FeSt 1-457 measured using the 30m IRAM tele- 
scope. Split asymmetric, double-peaked profile shapes characterize the CS and HCO + lines 
but not the C 18 line which is single-peaked indicating that the former emission lines are 
likely very optically thick and self-reversed. The sense of the asymmetry suggests expansion 
motions of the outer cloud layers. 
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Fig. 2. — Map of emission-line spectra from across FeSt 1-457. Solid lines are HCO + emission 
spectra and dotted lines are those of C 18 0. The map is centered at 0:2000 — 17 h 35 m 47.5 s , 
^2000 — -25°33'2.0". Double-peaked, self-reversed and asymmetric profiles are evident in 
the HCO + lines across the globule. In all positions the sense of the profile asymmetry is 
indicative of expansion motions. 
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Fig. 3.— Map of CS ( J = 2 -> 1) and C 18 ( J = 1 -> 0) emission- line spectra from 
FeSt 1-457. Solid lines are CS emission spectra and dotted lines are those of C 18 0. The 
map is centered at 0:2000 — 

17^35^47.5^ ,5 2000 = -25°33'2.0". Asymmetric, self-absorbed CS 
profiles indicative of expansion are evident across the map 
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Fig. 4.— Left. Peak velocity distribution for C 18 (J = 1 -> 0) emission 
Peak velocity distribution for C 18 ( J = 2 — > 1) emission spectra. 
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Fig. 5. — Line width distribution of C 18 ( J = 1 — > 0) spectra. The white cross denotes the 
position of the dust extinction peak. 
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Fig. 6. — Comparison of the spatial distributions of molecular emission shown as contours 
with a map of visual extinction constructed with the same angular resolution and shown as 
grey scale. Contour levels are given in J T* A dv with units of K km s _1 (C 18 0: 0.2 to 2.0 
by 0.2; N 2 H+: 0.1 to 0.8 by 0.1; HCO+: 0.2 to 1.2 by 0.2). The extinction scale is shown 
in the bar on the right of the figure and the extinction ranges from approximately 6-43 
magnitudes. For clarity the grey-scale contrast has been adjusted to emphasize the dust 
distribution in the center of the core. 
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Fig. 7. — {Left) C 18 J=l-0 integrated intensity as a function of visual extinction for the 
entire FeSt 1-457 dark cloud. {Right) N2H + J=l-0 integrated emission as a function of visual 
extinction. In all plots the data are presented as open squares with error bars while solid 
traces are the result of linear fits to the data. The extinction and molecular line observations 
in both plots were obtained at the same angular resolution. 
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Table 1. Line parameters. 



Line 


z/ a 


HPBW b 


Av c 




(GHz) 


(") 


(km s _1 ) 


CS (J = 2 -> 1) 


97.980968 


25 


0.030 


HCO+ (J= 1 -> 0) 


89.188512 


28 


0.033 


C 18 (J = 1 -> 0) 


109.78218 


22 


0.027 


C 18 (J = 2 -> 1) 


219.560319 


11 


0.013 


N 2 H+ ( J = 1 -> 0) 


93.173178 


26 


0.031 



a Line rest frequency 
b Half power beam width 
c Velocity resolution 
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Table 2. Results of velocity gradient fit. a 



Line 




v b 
(km s _1 ) 


dv/ds c 
(m s" 1 arcsec -1 ) 


dv/ds at 160 pc 
(km s _1 pc -1 ) 


9 d 

(°) 


C 18 ( J = 1 - 


,0) 


5.81 ±0.001 


0.73 ±0.012 


1.5 


231.5 ± 1.0 


C 18 ( J = 2 - 


-1) 


5.79 ±0.001 


0.70 ±0.022 


1.4 


238.9 ± 1.9 


N 2 H+ ( J = 1 - 


-0) 


5.79 ±0.003 


1.13 ±0.13 


2.3 


249.8 ± 7.5 



a Errors quoted are la uncertainty 
b Systemic velocity 
c Magnitude of velocity gradient 

d Direction of velocity gradient measured East of North 



